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We present molecular dynamics friction calculations for confined hydrocarbon solids with molec- 
ular lengths from 20 to 1400 carbon atoms. Two cases are considered: (a) polymer sliding against 
a hard substrate, and (b) polymer sliding on polymer. We discuss the velocity dependence of the 
frictional shear stress for both cases. In our simulations, the polymer films are very thin (~ 3 nm), 
and the solid walls are connected to a thermostat at a short distance from the polymer slab. Un- 
der these circumstances we find that frictional heating effects are not important, and the effective 
temperature in the polymer film is always close to the thermostat temperature. 

In the first setup (a), for hydrocarbons with molecular lengths from 60 to 1400 carbon atoms, the 
shear stresses are nearly independent of molecular length, but for the shortest hydrocarbon C20II42 
the frictional shear stress is lower. In all cases the frictional shear stress increases monotonically 
with the sliding velocity. 

For polymer sliding on polymer [case (b)] the friction is much larger, and the velocity dependence 
is more complex. For hydrocarbons with molecular lengths from 60 to 140 C-atoms, the number 
of monolayers of lubricant increases (abruptly) with increasing sliding velocity (from 6 to 7 layers), 
leading to a decrease of the friction. Before and after the layering transition, the frictional shear 
stresses are nearly proportional to the logarithm of sliding velocity. For the longest hydrocarbon 
(1400 C-atoms) the friction shows no dependence on the sliding velocity, and for the shortest hydro- 
carbon (20 C-atoms) the frictional shear stress increases nearly linearly with the sliding velocity. 
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1. Introduction 

Friction between solids is a very important phe- 
nomenon in biology and technology [l[ and it is very 
common in nature. Static friction always involves the 
coexistence of different metastable configurations at mi- 
croscopic level. When one surface slides on the other at 
low speed, first there is a loading phase during which the 
actual configuration stores elastic energy. Then, when 
the stored energy is large enough, an instability arises 
0,13, II: the system jumps abruptly to another configu- 
ration and releases elastic energy into irregular heat mo- 
tion. The exact way of how the energy is dissipated usu- 
ally does not influence the sliding friction force, provided 
that the dissipation is fast enough to happen before the 
next sliding event. 

There are many possible origins of elastic instabili- 
ties, e.g., they may involve individual molecules or, more 
likely, groups of molecules or "patches" at the inter- 
face, which have been named stress domains [EHSIll- 
Since the local rearrangements usually occur at different 
times in an incoherent manner, at the macroscopic scale 
the sliding motion may appear smooth without stick- 
slip oscillations. However, this is always a result of self- 
averaging, and at the atomistic level stick-slip motion will 



almost always occur (except for incommensurate systems 
with weak interactions). Moreover, at least at zero tem- 
perature, the friction force does not vanish in the limit 
of sliding speed v 0, but it tends to some finite value 
which depends on the average energy stored during the 
loading events and the atomic slip distance. 

A logarithmic velocity dependence of the frictional 
shear stress was observed in many experiments but usu- 
ally for sliding at low velocities (up to « 20 /im/s) 
Q. Still, in some experiments the logarithmic velocity 
dependence of the frictional shear stress was observed 
also for higher velocities (up to 1 mm/s) Com- 
paring experimental results with the existing theoretical 
models SL, JJJ shows that the logarithmic dependence 
could be well described in the models accounting for ther- 
mal activation effects. Thus, thermal fluctuations may 
induce jump of atoms (or rather group of atoms) at the 
sliding interface from one equilibrium position to the next 
one along the reaction path. The resulting stress-aided 
thermally activated effect leads to a logarithmic increase 
of friction with the velocity at low velocities. Thermal 
activation is more efflcient at low velocities, where the 
system spends long time in each potential well and, con- 
sequently, the probability to thermally activate the pro- 
cesses of atoms hopping is higher. 
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When a polymer (or a long-chain alkane) is sheared 
between two surfaces the shear stress often does not de- 
pend hnearly function of the logarithm of the sliding ve- 
locity This has been established both experimentally 
1^ . 13, 14 1 and in simulations jl5j. This observation may 
be due to the interdiffusion of chain segments between the 
polymer layers by the long molecules. Thus, the forma- 
tion of "bridges" is presumably affected by a change in 
the sliding velocity and in certain regimes it is possible to 
have decreasing shear stress with increasing sliding veloc- 
ity due to the disappearing of bridges across the sliding 
interface 14 , 1^ . This picture was first suggested in the 
context of rubber friction by Schallamach[245. 

In this paper we present molecular dynamics friction 
calculations for confined hydrocarbon solids with molec- 
ular lengths from 20 to 1400 carbon atoms. Two cases 
are considered: (a) polymer sliding against a hard sub- 
strate, and (b) polymer sliding on polymer. We discuss 
the velocity dependence of the frictional shear stress for 
both cases. We compare results obtained at room tem- 
perature and very low temperature (approaching K) . In 
the latter calculations no thermal activation can occur. 

In our simulations, the polymer films are very thin 
3 nm), and the solid walls are connected to a ther- 
mostat at a short distance from the polymer slab. Under 
these circumstances we find that frictional heating effects 
are not important, and the effective temperature in the 
polymer film is always close to the thermostat tempera- 
ture. In most practical situations the temperature is not 
fixed at planes close to the interface and during sliding at 
high enough velocities for a long enough time, the local 
temperature at the sliding interface may be so high as to 
locally melt the polymer surfaces. The physical processes 
occurring at the sliding interface in these cases may be a 
combination of the effects studied in this paper and the 
influence of the increased temperature. 



2. The model 

In this paper we present computer simulation results 
about the frictional behavior of linear hydrocarbons un- 
der applied pressure. Our model is similar to those de- 
scribed in Refs. SEHIli, 

but we review its main 
features here. We consider a block and a substrate with 
atomically flat surfaces separated by a polymer slab. Two 
cases are considered: (a) polymer sliding against a hard 
substrate which we will denote as "metal" for simplic- 
ity (the metal-polymer case), and (b) polymer sliding on 
polymer (the polymer-polymer case). 

The solid walls are treated as single layers of "atoms" 
bound to rigid flat surfaces by springs corresponding 
to the long-range elastic properties of 50 A thick solid 
slabs. For the case of sliding of polymer on "metal" , all 
molecules are adsorbed on the block surface only due to 
different parameters of interaction of alkane molecules 




FIG. 1: Snapshot picture of the C100H202 polymer slab at 
the sliding velocity n = 10 m/s and background temperature 
T — 300 K. (a) The molecules are (arbitrarily) colored in 
order to better observe the shear alignment of the chains, (b) 
The same as in (a) but with atoms presented as points in 
order to observe layering in the system. Seven monolayers of 
molecules are clearly seen. 



with the walls, whereas for the case of sliding of polymer 
on polymer about half of the molecules adsorbed on the 
block surface and half on the substrate surface. The block 
with adsorbed polymer slab was put into contact with the 
substrate surface in the first case and two solids with ad- 
sorbed polymer slabs were put into contact in the second 
case. When the temperature was equal to the thermo- 
stat temperature (usually 300 K) everywhere we started 
to move the upper block surface. We also conducted cal- 
culations for the temperature of solid walls equal to K 
in order to compare the sliding friction behavior at 300 
K with that in the absence of thermal activation at T = 
K. The temperature was also varied from 300 K to 550 
K to study the effect of melting on the shear stress. 

Linear alkanes C„Il2n+2 (with n ranging from 20 to 
1400) were used as "lubricant" in the present calcula- 
tions. The CH2/CII3 beads are treated in the united 



atom representation [20|, |21| . The Lennard- Jones poten- 



tial was used to model the interaction between beads of 
different chains 



Uir) = 4eo 



(?)"-"(?)* 



(1) 



and the same potential with modified parameters (ei, ri) 
was used for the interaction of each bead with the sub- 
strate and block atoms. The parameters were eo = 
5.12 meV for both the interior and the end beads, and 
To = 3.905 A and a = 1 in all cases. For the inter- 
actions within the Cnil2n+2 molecules we used the stan- 
dard OPLS model 2^, 21|, including flexible bonds, bond 
bending and torsion interaction, which results in bulk 
properties in good agreement with experimental data. 

For polymer sliding on polymer we need the polymer- 
metal bond to be so strong that no slip occurs at these in- 
terfaces. This is the case with ri = 3.28 A, ei = 40 meV 
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FIG. 2: The dependence of the shear stress on the shding 
velocity for the polymers from C20H42 to C1400H2802 sliding 
on "metal" at the normal pressure 10 MPa. The temperature 
300 K. The dotted line is a linear fit to the C20H42 curve. 



and a = 3. We also did some simulations with a — 2, but 
in this case some slip was observed at the polymer-metal 
interface. For sliding of polymer on "metal" we used the 
same parameters as above for the polymer-block inter- 
action (strong adsorbates interaction) but with a — I. 
For the polymer-substrate interaction we used a = I and 
£1 = 10 meV {weak adsorbate interaction). 

The choice of higher values of ei compared to eo re- 
flects the stronger (van der Waals) interaction between 
the beads and the "metal" surfaces than between the 
bead units of different lubricant molecules (this stronger 
interaction results from the higher electron density in the 
metals). The lattice spacings of the block and of the sub- 
strate are a = b = 2.6 A. 

We used linear alkane molecules with the number of 
carbon atoms 20, 60, 100, 140 and 1400 as lubricant. 
The number of C100H202 molecules was equal to 200. The 
number of C20H42 molecules was equal to 1000. This gave 
from 6 to 8 monolayers of lubricant molecules between 
the solid surfaces. The (nominal) squeezing pressure pq 
was usually 10 MPa. 

As an illustration, in Fig. [T] we show the contact be- 
tween a flat elastic block (top) and a flat elastic substrate 
(bottom). The polymer slab (^ 30 A thick) is in between 
them. Only the interfacial block and substrate atoms and 
polymer atoms are shown. 



3. Polymer on "metal" 

In Fig. 12] we show the dependence of the shear stress 
on the sliding velocity for the polymer slabs sliding on 
"metal" at the appHed pressure p — 10 MPa. For C20H42 
the lubricant behaves liquid-like with the shear stress 
nearly proportional to the sliding velocity (see dashed 
curve in Fig. In particular, as the sliding velocity ap- 
proaches zero the shear stress for C20H42 becomes very 
small, reflecting the fact that this polymer is close to the 
liquid state: the C20H42 (eicosane) system is at 27°C 
(the temperature of the thermostat), which is very close 
to the melting point (« 37°C) of this polymer. We also 
performed simulations at the (thermostat) temperatures 
200 K and K, and in these cases the friction is consid- 
erably higher (~ 9.6 times higher at the temperature K 
for the sliding velocity v = 0.3 m/s). 

For longer-chain hydrocarbon molecules we found that 
the polymer films behave solid-like at 300 K, with non- 
vanishing kinetic friction as the sliding velocity ap- 
proaches zero, and that the shear stress depends non- 
linearly on the sliding velocity. This is due to the higher 
melting point of the longer-chain polymer systems. For 
the C100H202 system the frictional shear stress is more 
than twice higher at T = K than for T = 300 K for the 
sliding velocity v — 0.3 m/s. This can be attributed to 
the absence of thermal activation at T = K. Atomistic 
stick-slip events occur at the sliding interface between 
the polymer and the substrate. The period of oscilla- 
tions of the shear stress as a function of x-coordinate of 
the block is exactly 2.6 A which is one lattice unit (of the 
substrate) in the sliding direction. At low sliding veloc- 
ities the shear stress increases greatly (much more than 
for the higher sliding velocities) when the temperature 
is decreased from 300 K to K. This is due to the fact 
that the probability to activate the processes of atoms 
hopping is higher at low velocities due to the longer time 
the system spends in each local potential well along the 
reaction coordinate. 

We emphasize the importance of the temperature (or 
thermal fluctuations) on the process of "going over the 
barrier" . Thus at zero temperature, the external ap- 
plied tangential force (or stress) alone pulls the system 
over the lateral pinning barriers, and this happens ev- 
erywhere simultaneously. At high sliding velocities ther- 
mal effect should be rather unimportant. However, for 
small sliding velocities, thermal fluctuations will be very 
important. In this case slip will not occur everywhere 
simultaneously, but small nanometer-sized interfacial re- 
gions of linear size D will be individually pinned and 
perform stress-aided thermally induced jump from one 
pinned state to another (local interfacial rearrangement 
processes). (Note that thermal effects can only become 
important for small (nanometer-sized D) regions, since 
simultaneous going-over-the-barrier everywhere requires 
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FIG. 3: The shear stress for Ci4oH[282 polymer sliding on 
metal at the sliding velocity v = 10 m/s and the normal 
pressure p = 10 MPa when ramping (step-like changing) the 
temperature T from 300 K to 550 K (heating) and from 550 K 
to 300 K (cooling). 



infinitely large energy for an infinite system, except, per- 
haps for an incommensurate interface.) This pro cess has 
been studied in detail both theoretically 22, 2^ 24, 25| 
and experimentally [sl. [26l. [27t . 

We now study the frictional behavior of the hydrocar- 
bon films C100H202 and C140H282, when the temperature 
is increased above the melting temperatures. In Fig. [3] we 
show the shear stress for the C140H282 polymer film (at 
the sliding velocity u = 10 m/s and the normal pressure 
p = 10 MPa) when increasing the thermostat tempera- 
ture T by steps equal to 50°C from 27°C to 277°C, and 
then decreasing it back to 27°C. During heating the shear 
stress decreases abruptly when the temperature is raised 
above the melting temperature Tm = 110°C for C140H282 
polymer. The frictional shear stress for the film in the 
liquid-like state is ^ 3 times lower than for the solid film 
just below melting. 

For temperatures well above the melting point the 
molecules in the center of the polymer film are disor- 
dered as expected for the liquid state of the lubricant. 
In Fig. |4] we show the density distribution, and the aver- 
age velocity Vx of the C-atoms along the distance between 
the substrate and the block (z-direction) for the C100H202 
polymer film. The results are for T = 27°C and 177°C. 
At T = 27° C the film is seven monolayers thick, but at 
T = 177°C, due to the thermal expansion, the system 
has eight monolayers. In the center the molecules are 
disordered and the layers of molecules blurred. 

For 27°C the average velocity changes abruptly at the 
substrate-lubricant interface (see Fig. 2^), so the slip oc- 
curs between the substrate and the first monolayer of 



molecules of the polymer film. Thus, the whole poly- 
mer film is bound to the block and moves with the av- 
erage velocity ~ 10 m/s, i.e. with velocity of the 
block. For T = 177°C most of the slip also occurs at 
the polymer-substrate interface (see Fig. [3)3), but a small 
slip (slip velocity w « 1 m/s) also occurs at the polymer- 
block interface. Thus, all monolayers of the polymer film 
move with the average velocity Vx ^ Q m/s. The slip at 
the polymer-block interface is due to the applied shear 
stress and thermal fluctuations. At the lower tempera- 
ture 27° C the thermal fluctuations are not strong enough 
to overcome the relatively large atomic corrugation at the 
polymer-block interface, which results from the relatively 
large ei = 40 meV for the interaction of each polymer 
bead units with the block atoms (strong adsorbate in- 
teraction). Friction is a stress-aided thermally activated 
process, and in the present case when the temperature in- 
creases from T = 27°C to 177°C, the shear stress drops 
by a factor of ^ 2. 

The effective corrugation of the interaction potential 
experienced by the molecules at the sliding interface is 
the most important parameter influencing the magni- 
tude of the friction and the dependence on the exter- 
nal (squeezing) pressure. Indeed, the fact that the lat- 
tice constant of the substrate is much smaller than the 
size of the polymer molecules (and the polymer persis- 
tence length), and also very different from the natural 
separation between the polymer molecules, implies that 
the effective corrugation of the interaction potential be- 
tween the polymer and the "metal" substrate will be very 
small, and this explains the small friction observed in this 
case, compared to the case when the slip occurs at the 
polymer-polymer interface. 



4. Polymer on polymer 

When a polymeric fllm is strongly attached to the block 
and substrate surfaces, sliding of the block will induce a 
shearing of the polymer fllm. For this case we have in- 
vestigated the influence of the sliding speed on the shear 
stress. Figure [5] shows the results obtained for a number 
of linear alkanes with chain length from 20 C-atoms to 
1400 C-atoms. Note that for C1400H2802 the shear stress 
is independent of the sliding velocity, whereas the C20H42 
is more liquid-like, with a shear stress approximately lin- 
early related to the sliding velocity. 

For the mid-sized molecules (with 60 to 140 C-atoms in 
the chain), for small and large velocities the shear stress 
is nearly proportional to the logarithm of the sliding ve- 
locity, with a small slope at low velocities and a larger 
one at large velocities. In the range of from 20 m/s to 
40 m/s the shear stress decreases with increasing sliding 
velocity. This is due to a layering transition, where the 
number of layers increases from 6 to 7. This is proved in 
Fig. [5] which shows the positions in the the z-direction 
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FIG. 4: The distributions of the average number and the 
average velocity (in the shding direction) of lubricant C 
atoms along the distance between the substrate and the block 
for C100H202 metal-polymer sliding at velocity of the block 
V — 10 m/s for the normal pressure p — 10 MPa for the 
temperature a) 27° C and b) 177° C. In the latter case the 
slip also occurs between the block and the last monolayer 
of lubricant molecules due to thermal fluctuations. The very 
left and very right maximum of the both density distributions 
correspond to the substrate and the block atom layers. 



of the layers for the C100H202 system as a function of the 
sliding distance d. At d = the sliding velocity of the 
block is changed from 100 m/s to 10 m/s and after some 
relaxation time period (which corresponds to the sliding 
distance d w 1000 A) the system abruptly switches from 
7 to 6 layers. At the same time the shear stress increases 
abruptly (see Fig. [B]). The latter is, at least in part, 
due to the decreased number of slip planes. The layering 
transition in Fig. [6] is reversible: increasing the velocity 
back to 100 m/s results in a return to 7 layers. 

As pointed out above, for the mid-sized molecules, for 
small and large velocities the shear stress is nearly pro- 
portional to the logarithm of the sliding velocity, with 
a small slope at low velocities and a larger one at large 
velocities. The logarithmic velocity dependence is ex- 
pected for thermally activated stress induced processes. 




sliding velocity (m/s) 



FIG. 5: The shear stress as a function of the sliding veloc- 
ity for all the investigated systems. The normal pressure is 
10 MPa. 



which predict the frictional shear stress 
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where k^T is the thermal energy (T is the temperature), 
Eb is an energy barrier (activation energy for some rear- 
rangement process involved in lateral slip, e.g., removal 
of polymer bridge) and vq is a reference velocity. Thus, 
the larger slope of the Uf (logw) relation for high slip veloc- 
ities, as compared to low slip velocities, can be explained 
by assuming that the energy barrier Eb is smaller in the 
more open structure, which prevail after the transition 
from 6 to 7 layers with increasing velocity. 

The change in the number of layers in the film is ther- 
mally activated, since no change in the number of layers 
(on the time scale of our simulations) occurs when the 
thermostat is at K. Experimental data suggest that 
the layering transition happens at lower sliding velocities 
when the normal pressure is decreased [l3j . This is in 
good accordance with our results as a decrease in pres- 
sure from 10 MPa to 3 MPa shifted the transition velocity 
of the C60H122 system from about 30 m/s to 20 m/s. 

We will now study the nature of the layering transi- 
tion in greater detail. In Fig. [7] we show the separation 
between the surfaces for each polymer system we have 
studied, as a function of the sliding speed. Note that the 
surface separation in the case of C20H42 and C1400H2802 
increases with increasing sliding velocity with a constant 
slope, but the number of polymer layers is constant (at 
seven) in these cases. For the mid-sized molecules there 
is an abrupt increase in the surface separation in the 
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FIG. 6: The shear stress and the positions in the the z- 
direction of the layers for the C100H202 system as a function 
of the shding distance. The sliding velocity of the block is 
changed from 100 m/s to 10 m/s at A. 



transition regime of figure [S) The number of layers is 
six before the transition (below 10 m/s) and seven after 
(above 40 m/s). 

Let us now study the number of molecular bridges be- 
tween the layers, which may be affected by the layering 
transition [ij, [3, [3, [3] . There is no exact definition of 
a molecular bridge, but we have chosen to define a bridg- 
ing atom by the fact that it does not belong to the same 
layer as the preceding atom in the molecule, and at least 
three of the preceding and three of the following atoms 
in the molecule are in two different layers, see Fig. [51 
The bridging phenomenon can be observed in the snap- 
shot shown in figure [51 The number of bridging atoms 
as a function of the sliding velocity is shown in figure [TOl 
Note that the C1400H2802 system has the same number 
of bridges through the whole velocity range. The C20H42 
system shows a linear dependence, whereas the mid-sized 
molecules seem to have a constant number before the 
transition regime, increasing to a higher level after this 
regime. The latter can be explained by the more open 
structure of the 7-layer systems, which makes it easier for 
polymer molecules from one layer to have segments ex- 
tending into other layers. Figure [H also shows that some 
atoms are outside the center line of the polymer layers 
without linking these. We define such a non-layer atom 
as one having a distance of at least 1.5 A to the center 
line of any layer (see Fig. [8]). We present the number of 
non-layer atoms as a function of the sliding velocity in 
figure [TTl 

Figure [TT] shows that the C1400H2802 system has an in- 
creasing number of non-layer atoms with increasing slid- 
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FIG. 7: The distance between the surfaces as a function of 
the sliding velocity for all the investigated systems. There is 
also an indication of the number of layers present at differ- 
ent surface separations. It was chosen to use a linear x axis 
to show the transition more clearly. The normal pressure is 
10 MPa. 
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FIG. 8: The definitions of bridge and non-layer atoms used 
in the present paper. 
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FIG. 9: Snapshot picture of a C100H202 system where ten 
molecules have been chosen to be shown at random. The rest 
of the atoms are reduced to points. The picture shows that 
the molecules have atoms in several layers. Some of these 
segments form bridges but others are just present in different 
layers without linking these. The normal pressure is 10 MPa 
and the sliding velocity is 0.3 m/s. 
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FIG. 10: The number of bridging atoms as a function of 
the sliding velocity for the investigated systems. A bridging 
atom is defined in the text and in figure [S] The dotted line is 
a linear fit to C20H42. The normal pressure is 10 MPa. 
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FIG. 11: The number of non-layer atoms as a function of 
the sliding velocity for the investigated systems. A non-layer 
atom is defined in the text and in figure |8] The dotted line is 
a linear fit to C2oIl42. The normal pressure is 10 MPa. 



ing velocity. In the case of C20H42 the number of non- 
layer atoms is proportional to the sliding velocity. The 
increase in the number of non-layer atoms is associated 
with the increase in the separation between the layers 
with increasing sliding velocity (see Fig. [7]), which al- 
low polymer segments in a layer to displace away from 
the layer-plane by a considerable distance without expe- 
rience a large repulsion from the nearby polymer layer. 
One may alternatively interpret the increase in the layer 
separation as resulting from the increased repulsion from 
the non-layer atoms with increasing sliding velocity. In 
this picture the increased number of non-layer atoms re- 
sults from the increased momentum transfer in collisions 
between atoms in two nearby layers as the sliding ve- 
locity increases: these collisions kick polymer segments 
away from the layer-plane. 

The C60H122, C100H202 and the C140H282 systems ex- 
hibit a slow increase in the number of non-layer atoms 
before the transition, and a faster one when the polymer 
film has increased to seven layers. Note that the C100H202 
system exhibits a maximum in the number of non-layer 
atoms around 10 m/s. This is caused by a strong fluc- 
tuation which sometimes was observed also for the 60 
C-atom and 140 C-atom systems (not shown). This is 
not unexpected as strong fluctuations often occur close 
to phase transition points (in this case a layering transi- 
tion). 

Let us now compare the transition observed in figure [5] 
with what happens during melting of the polymer film. 
We have investigated the behavior of the film when it is 



melted by raising the temperature from 300 K to 450 K; 
see figure [l2l The figure shows the shear stress and the 
distance between the surfaces as a function of the sliding 
time. For both 3 and 10 m/s the shear stress decreases 
to about 10 MPa, a much lower level than observed in 
figure [5l At 3 m/s the film keeps its 6 layers as the 
melting increases the distance between the surfaces to 
about 27.4 A, a result that is consistent with figure [T] 
The combination of a higher sliding velocity and melting 
increases the distance between the surfaces to about 29 A 
when the sliding velocity is 10 m/s. The film then passes 
to 7 layers as predicted by figure [71 From this we deduce 
that the layering transition induced by the increase in 
the sliding velocity is not associated with a melting of the 
film. The transition seems to occur at a certain surface 
separation and the shear stress is not affected by it in 
the melted state whereas this is the case when the sliding 
velocity is increased. 

Let us study the velocity profiles of the systems. Figure 
fT3l shows the cumulative velocity probability distribution. 
Note that for the C1400H2802 system nearly all the slip 
occurs at one interface at the center of the polymer film. 
Thus, this system can be considered as two polymer slabs 
sliding against each other. This picture also explains the 
independence of the shear stress for C1400H2802 to the 
number of bridges (figure fTO]) and non-layer atoms (figure 
fTTj) . since these are presumably mostly appearing inside 
the polymer slabs, and thus do not influence the interfa- 
cial slip. The C20H22 molecules have their sliding planes 
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FIG. 12: The shear stress and the distance between the sur- 
faces as a function of the sliding time when the temperature 
is increased from 300 K to 450 K. The melting point is 390 K. 
The sliding velocities are 3 and 10 m/s. The normal pressure 
is 10 MPa. During melting the film increases its number of 
layers from 6 to 7 at 10 m/s whereas it stays at 6 layers at 3 
m/s. The system is C100H202. 

distributed over the whole thickness of the film. This 
is expected for a liquid-like flow, but figures [TO] and [11] 
indicate that this also could be a result of an increasing 
number of bridges and/or non-layer atoms. The mid- 
sized molecules in figure [T3] have a velocity profile in be- 
tween the longest molecules and the shortest ones, and 
can be considered as a so-called plug flow, where the out- 
ermost layers are pinned to the surfaces of the block and 
the substrate. 

In Fig. [TU at d = A (where d is the slip distance) 
the sliding velocity of the block is changed from 100 m/s 
to 10 m/s. Note that as the fllm relaxes the number of 
bridges is quite constant whereas the number of non-layer 
atoms is decreasing slightly and so is the shear stress. 
When the transition regime is reached the number of 
non-layer atoms suddenly increases as the central lay- 
ers become mobile and (in the z-direction) diffuse. The 
increase in shear stress follows this increase in the num- 
ber of non-layer atoms whereas the number of bridges 
decreases to a lower constant level as the number of lay- 
ers goes from seven to six. This can be understood as the 
density of molecules in the layers is smaller for the 7-layer 
state than for the 6-layer state, and hence the ability for 
chain molecules to (due to a fluctuation) interdiffuse and 
form bridges will be largest in the 7-layer state. 

The shift from six to seven layers increases the num- 
ber of slip planes by one. At the same time the density 
of bead-units in the layers decreases so that the space 
between the bead units in each layer increases. The in- 
creased space reduces the energetic barriers for polymer 




-2 2 4 6 8 10 12 
sliding velocity (m/s) 



FIG. 13: Shding velocity of atoms as a function of the number 
of atoms sliding at less than this specific sliding velocity. The 
number of atoms is accumulated so that the total reach 20000 
for C20II42 and CiooIl202 whereas the C1400II2802 system only 
has 19600 atoms. The sliding velocity of the block is 10 m/s 
and the normal pressure is 10 MPa. 

segment rearrangement processes, and results in the for- 
mation of more cross-links (see figure [TT|) . At low veloci- 
ties the shear stress increases slowly with the increase of 
the sliding velocity, but after the transition the slope be- 
comes steeper (see figure[S]). This can again be attributed 
to reduced energetic barriers for polymer rearrangement 
processes. 

5. Summary and conclusions 

We have presented results of molecular dynamics cal- 
culations of friction for two solids separated by a w 3 nm 
thick polymer film. Two types of systems were consid- 
ered: (a) a polymer film pinned to one of the solid sur- 
faces and sliding at the other solid surface (the "metal" - 
polymer case). The second case (b) was with the polymer 
layers pinned to both solid surfaces and shearing at the 
polymer-polymer interface (the polymer-polymer case). 
We used linear alkane molecules with the number of car- 
bon atoms 20, 60, 100, 140 and 1400. 

The frictional shear stress for the polymer-polymer 
systems is much higher than for the "metal" -polymer 
systems. This is due to the same size of the atoms or 
molecules on both sides of the slip-plane for the polymer- 
polymer case, resulting in strong interlocking (as for a 
commensurate interface), while the "metal" -polymer in- 
terfaces are incommensurate (the lattice constant of the 
"metal" substrate is different from the distance between 
atoms of the lubricant molecules) . 
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FIG. 14: The number of bridging and non-layer atoms and 
the shear stress as a function of the sliding distance for 
CiooH202. The sliding velocity of the block is changed from 
100 m/s to 10 m/s at A. 



We have studied the velocity dependence of the fric- 
tional shear stress for both cases. In the first setup 
the shear stresses are relatively independent of molecular 
length. For the shortest hydrocarbon C20H42 the fric- 
tional shear stress is lower and increases approximately 
linearly with the velocity. 

For polymer sliding on polymer the friction is signifi- 
cantly larger, and the velocity dependence is more com- 
plex. For the longest molecules (1400 carbon atoms) the 
shear stress is independent of the sliding velocity as the 
sliding occurs primarily at one interfacial slip plane. The 
shortest molecules again exhibit liquid-like sliding with 
the shear stress being approximately proportional to the 
sliding velocity. The mid-sized molecules (60 to 140 C- 
atoms) show a slightly increasing shear stress at low ve- 
locities, and a faster increase at high sliding velocities. 
Between these regimes there is a transition with a de- 
crease in the shear stress with increasing sliding velocity. 

The mechanism behind this behavior seams to be a ki- 
netic phase transition involving a change in the number of 
layers in the film, which introduce new slip planes. This 
decreases the shear stress abruptly. Further increase of 
the sliding velocity will increase the shear stress rapidly, 
which we attribute to the interaction between the layers 
via non-layer atoms. 
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